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ABSTRACT 

O ■ 

>/") . We report on the changes of the Sun's subsurface stratification inferred from 

helioseismology data. Using SOHO/MDI (SOlar and Heliospheric Observa- 
tory/Michelson Doppler Imager) data for the last 9 years and, more precisely, 
q the temporal variation of f-mode frequencies, we have computed the variation 

& ', of the radius of subsurface layers of the Sun by applying helioseismic inversions. 

We have found a variability of the "helioseismic" radius in antiphase with the 
j> ! solar activity, with the strongest variations of the stratification being just below 

^ , the surface around O.995i? . Besides, the radius of the deeper layers of the Sun, 

between O.975i? and O.99i? changes in phase with the 11-year cycle. 

Subject headings: Sun: helioseismology — Sun: oscillations — Sun: activity — 
Sun: interior 



1. Introduction 

For the last decades, temporal variations in the solar radius has been a controversial 
subject. Indeed, measurements made with the solar astrolabe by F. Laclare in France (La- 
clare et al. 1996) and the Brazilian team (Reis Neto et al. 2003) showed a variation of the 
solar radius in antiphase with the solar activity cycle, while F. Noel in Chile (Noel 2004) 
using a similar instrument reported a variation in phase with the solar activity. By using 
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solar continuum intensity images obtained with the Michelson Doppler Imager (MDI) on 
board SOHO, Emilio et al. (2000) found a variation of the solar radius of 8.1 ± 0.9 mas/yr, 
but in a recent and more complete study of these images, Kuhn et al. (2004) reported no 
evidence of solar-cycle visible radius variations between 1996 and 2004 at any level above 
7 milliarcseconds: this is significantly lower than any variation reported from ground-based 
observations. Other measurements made by the Solar Disk Sextant (SDS) experiment on- 
board strastospheric ballons (Sofia et al. 1994; Thuillier et al. 2005) reported a variation of 
the solar radius in antiphase with the solar radius. The solar radius can be also determined 
by helioseismic methods (Schou et al. 1997). This, so-called "seismic" radius, is related to 
the subsurface density stratification and can be compared with the "photospheric" radius, as 
inferred by astrolabe ground-based measurements, for example, only by using solar models. 
The "seismic" radius probes the sub-photospheric layers up to a depth of about 15 Mm. 
So the use of solar f-mode frequencies to infer the seismic radius is important for searching 
for physical changes occuring beneath the photosphere. Results are, however, conflicting. 
Dziembowski et al. (2001) reported rates of the seismic radius change ranging from -3 to 1 
km/year during the rise of cycle 23. Recently, Dziembowski & Goode (2005) reexamined the 
issue by using SOHO/MDI f modes for 1996-2004, and were unable to obtain stable results 
for the seismic radius, and thus they concluded that the observed variations of the f-mode 
frequencies should be explained by the direct effect of magnetic field. However this led to 
unrealistically strong random field hidden below the surface. So, the question whether or 
not there is a radius variation with the activity cycle is still debated. Moreover if there are 
such variations at the surface, where is their origin? In this letter, we consider the issue 
of solar radius and determine solar radius variations with time for various layers below the 
surface. The main difference from the previous investigations is that we do not assume an 
uniform change of the radius of subsurface layers, but allow variations of displacement of 
these layers with depth. We show that a stable solution does exist in this case, and find 
evidence for temporal solar radius variations in the sub-photospheric layers above O.97_R 
with an oscillation, in antiphase with the solar cycle above O.99-R and in phase below. Our 
results show the localization of these variations in the upper convective zone. If we extrapo- 
late these results up to the surface, we find a radius change of about 2 km in antiphase with 
the solar cycle: here, we have to keep in mind that without high-1, we cannot constrain the 
surface radius better, and that in reality, this variation at the surface can be larger provided 
it is more localized. 
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2. Data 

We used frequencies of solar oscillation modes from 72-day MDI observing runs, com- 
puted by J. Schou 1 . We selected only f modes for the observing period 1996-2004. Each file 
has a different number of f modes, so we extracted only common modes from these files and 
obtained a total of 151 modes ranging from the angular degree I = 125 to I = 285. We have 
computed the relative frequency changes bvjv for each degree and calculated the average 
over each year and binned every 20 /zHz. The reference year has been chosen to be 1996, 
near the minimum of the solar cycle 2 . The averaging over 1-year of the data allowed us to 
avoid the instrumental 1-year variations. The results are compiled on Fig. 1. The errorbars 
are not plotted for clarity of the graph. 

The different curves plotted on this figure shows an evidence of variations in the f- 
mode frequencies over the solar cycle already pointed out by Dziembowski et al. (2001). We 
emphasize an intriguing phenomenon at higher /, above v = 1600 /xHz, where the frequency 
difference can change of sign, becoming then negative, in the descending part of the solar 
cycle. It is particularly puzzling that the sharp frequency decrease above 1600 yuHz continues 
in the declining phase of the solar cycle whereas the frequencies of the lower-frequency modes 
return to their previous solar minimum values. We don't know how to interpret this change 
but assume that it may come from a variation in the near-surface turbulence, which can 
affect the frequencies of the f modes confined just below the surface, in a zone close to the 
leptocline (from the greek /eptos=fine), thin transition layer between the upper convective 
zone and the photosphere (Godier & Rozelot 2001; Rozelot & Lefebvre 2003). We used these 
computed frequency differences (without the binning over v) between / = 150 and 250, which 
are measured most reliably, as input parameters in our helioseismic inversion presented in 
the next section. 



3. Helioseismic inversion of f modes to infer solar radius variations 

3.1. Mathematical formalism 

As a starting point for our inversion, we used the equation derived by Dziembowski 
& Goode (2004) who established a relation between the relative frequency variations 5v/v 
for f-mode frequencies and the associated Lagrangian perturbation of the radius 5r/r of 



1 http://quake. stanford.edu/~schou/anavw72z/ 
2 For this year, data begin from May I s *. 
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subsurface layers: 

'5v\ 31 f g Sr 



, dl-— (1) 

v J x 2cu 2 I J r r 

where / is the degree of the f modes, / is the moment of inertia, oo the eigenfrequency and g 
the gravity acceleration. 

This equation leads to the asymptotic relation used in the previous determinations of 
the solar radius using f modes ^ = — (Schou et al. 1997), assuming that y is constant 
with depth. However, we don't make this assumption, and infer y as a function of r. Eq. 1 
can be rewritten as 

- = / K x -dr 2 
v , . n r 



where the kernel Ki is expressed as 



Ki = ^ 

1 O, ,2 



3/ 



gr (3) 



p being the density and the mode eigenfunction. 



To compute each kernel Ki, we used the model S (Christensen-Dalsgaard et al. 1996) 
calibrated to the seismic radius of Schou et al. (1997) (R Q = 6.9568 x 10 5 km). Fig. 2 shows 
three examples of the kernels at / = 150, 200 and 250. The method used to invert Eq. 2 is 
the standard Regularised Least-Square technique (Tikhonov & Arsenin 1977), appropriate 
for this kind of ill-posed problems. In this method, one has to minimize the following relation 

E = — \ Ki y dr -fi) +a y 2dr 

j a i \Jo J Jo 







where y — y, fi — (y) r <J\ the relative incertainty for each fi, a and (3 the regular isat ion 
parameters. 



3.2. Results 

In the inversion process, we restricted the data set using only modes with / ranging from 
150 to 250 measured more accurately than the other modes. For each year, the minimization 
of Eq. 4 leads to a solution plotted on Fig. 3. This solution integrated through Eq. 2 yields 
the relative frequencies variations versus the degree / in comparison with the real data. The 
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reconstructed frequencies fit pretty well the trend of the real relative frequencies within the 
errorbars (Fig. 4). Fig. 5 shows the AR estimated at the surface using R Q = 6.9568 x 10 5 
km as a reference in the model. Fig. 6 shows the averaging kernels (Ory & Pratt 1995) 
which illustrate the resolving power of these inversions. It is quite clear that the spatial 
resolution (FWHM) is about O.OO5i? in the region between O.98i? and O.998i? . Like in 
other inverse problems, this means that the variations of the smaller scale and outside this 
region cannot be resolved without additional constraints. In Fig. 7 we give three examples of 
test inversions with artificial f-mode frequency data (applying the observed error estimates), 
which illustrate the accuracy and limitations of these inversions. 

The main characteristics of our solution are: 

1. Fig. 3 shows no significant changes in the variation of the subsurface layers depth 
below 0.97i? o . In this layers, our inversions loose the spatial resolution as follows from 
Fig. 6. So, with the currently available data, it is not possible to measure localized 
variations below O.97i? . However, if there were an uniform change of the radius in 
subsurface layers (including O.97i? and below) of the order of 1-3 km as discussed by 
Dziembowski et al. (2001), it could be detected by the inversion procedure as illustrated 
in the top panel of Fig. 7 (no assumption was made about the functional form of the 
solution). 

2. Fig. 3 shows non-monotonic changes in the stratification with the inner layer (below 
O.99i? ) moving up during the increase of activity (compression) and the outer layer 
(above O.99i? ) moving down (relaxation). The precise localization of these layers is 
uncertain because have a characteristic width of about O.OO5i? . The test inversion in 
Fig. 7 (middle panel) shows that this uncertainty can be about the half-width of the 
averaging kernels, about O.OO3i? . 

3. Fig. 5 estimates that the near-surface variations are in antiphase with the solar cycle 
with an amplitude of the order of 2 km. However, the sensitivity of our inversions 
is quite low at the surface, and localized variations of the surface radius may not be 
detected (see a test inversion in bottom panel of Fig. 7). High-degree f-mode data are 
required to improved the surface estimates. 

Note that if we used all the modes available (i.e. 151 modes), we cannot find a stable 
solution that fits the last part of the curves above v = 1600 fiHz (see Fig. 1). We suppose 
that the behavior of the curves in this range of frequencies could be due to the influence of 
turbulence and magnetic fields very near the surface. 
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4. Discussion 

First of all, Fig. 3 shows a temporal variation of the solar radius in the subsurface layers 
during the solar cycle. This oscillation is centered around r = O.99i? with a width of about 
O.O3i? . This anomaly is composed of two parts: the first, below r = O.99_R , is in phase 
with the solar cycle and has a maximal amplitude of about 10 km; the second part above 
r = O.99i? , is in antiphase with the solar cycle and reaches larger amplitudes, of about 26 
km. These variations indicate the presence of a changing with the solar cycle of a physical 
structure that could be described as a very thin transition layer, siege of the variation of the 
solar radius. This transition layer is located here approximately at O.99i? and it is linked 
to changes in the upper convective zone caused by magnetic fields. We estimate a variation 
of the seismic radius at the surface of about 2 km at the maximum of the cycle 23. However, 
the surface radius is poorly constraint with the currently available set of medium-1 f-mode 
frequencies. 

Our helioseismic inversion computations have put in evidence a variability of the solar 
radius in the subsurface layers which are extended to the surface. The results are gener- 
ally consistent with previous conclusions that solar-cycle variations in the solar radius are 
confined to the outermost layers of the Sun (Antia & Basu 2004; Dziembowski & Goode 
2005). This variability is localized in a double-structure layer centered at O.99i? : in the 
deeper part, between O.97i? and O.99i? , the radius varies in phase with the solar cycle, 
whereas this is opposite in the upper part above O.99i? , where the variability become in 
antiphase. Thus we confirm the fact that the surface layers of the Sun are shrinking during 
the ascending phase of the solar cycle and is relaxing after the maximum. However, these 
changes are not uniform with depth. In a near future, it would be interesting to inspect 
more in details the behavior very close to the surface by looking at changes of the higher 
degree modes, above / = 250, where a second thin layer may take place. 

As a conclusion, we would like to emphasize our most significant result: the change in 
radius goes from being in phase with the solar cycle in the deeper layers to out of phase 
in the shallower layers with a transition at O.99i? . This result could eventually lead to a 
deeper understanding of the physics behind the changes. 

This work utilizes data from SOHO/MDI and we thank J. Schou for providing the 
frequencies files. 
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Fig. 1. — Average relative frequency differences in f-mode (Su/u) as a function of (v>), average 
frequencies binned every 20 /iHz. The reference year is 1996 and the errorbars have not be 
plotted for clarity of the graph. The f modes chosen for our study have frequencies between 
the limits represented by vertical lines. 




Fig. 2.— Example of kernels K t at I = 150, 200 and 250. 
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Fig. 3. — (Sr) as a function of the fractional radius x = r/R Q , obtained as a solution of the 
minimization of Eq. 4. Notice the behavior of the curve near x = 0.99. The errors are the 
standard deviation after average over a set of random noise added to the relative frequencies 
reconstructed in Fig. 4. 



- 11 - 




Fig. 4. — For each year, (Sv/u) as a function of the degree /. The reference year is 1996 and 
the errors are the relative incertainties. The solid curve is the results of direct integration of 
Eq. 2 providing the solution of Fig. 3. 
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Fig. 5. — Top: temporal variation of AR near the solar surface at r = R Q from the solution 
of Fig. 3. Bottom: Variation of the sunspot number for the same period. The variation of 
the seismic radius at the surface is found in antiphase with the solar cycle with an amplitude 
of about 2 km. It is important to keep in mind that without high-1, we cannot constrain the 
surface radius better, and that in reality, this variation at the surface can be larger provided 
it is more localized. 
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Fig. 6. — Amplitude of the averaging kernels versus the fractional radius 



-14- 




9.5 
40 

1 20 

to 



40 

I 20 

to 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i- 
0.965 0.97 0.975 0.98 0.985 0.99 0.995 1 



T 



T 



TT 



T 



_L 



_L 



. ■■ * ■* I ^ 



J_L 



_L 




0.965 0.97 0.975 0.98 0.985 0.99 0.995 1 



T 



T 



TT 



T 



T 



J_L 



m-t-l-m-iwrf 




0.965 0.97 0.975 0.98 0.985 0.99 0.995 1 

x=r/Ro 

Fig. 7. — Inversion tests with artificial data; from top to bottom: the initial 5r/r is a 
constant, a gaussian with a width of 0.005 and centered on x=0.99, a gaussian with a width 
of 0.005 and centered on x=l. In all cases, the solid line is the initial data and the dashed 
curve is the result after inversion. 



